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r^ ' Abstract. We review recent developments in the study of excited baryon spectroscopy 

in lattice QCD. After introducing the basic methods used to extract masses from cor- 
relation functions, we discuss various interpolating fields and lattice actions commonly 
used in the literature. We present a survey of results of recent calculations of excited 
baryons in quenched QCD, and outline possible future directions in the study of baryon 
spectra. 
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en ■ 1 Introduction and Motivation 

^^ I One of the primary tools used for studying the forces which confine quarks inside 

^-r ■ hadrons, and determining the relevant effective degrees of freedom in strongly 

^s I coupled QCD, has been baryon and meson spectroscopy. This is a driving force 

behind the current experimental N* programme at Jefferson Lab, which is ac- 
cumulating data of unprecedented quality and quantity on various N — > N* 
transitions [1]. The prospects of studying mesonic spectra, and in particular the 
O , role played by gluonic excitations, has been a major motivation for future facil- 

2 ■ ities such as CLEO-c [2] , the anti-proton facility at GSI (PANDA) [3] , and the 

Hall D programme at a 12 GeV energy upgraded CEBAF [4]. With the increased 

precision of the new N* data comes a growing need to understand the observed 

spectrum within QCD. QCD-inspired phenomenological models, whilst success- 

5^ ' ful in describing many features of the N* spectrum [5], leave many questions 

^ ! unanswered. 

One of the long-standing puzzles in baryon spectroscopy has been the low 
mass of the first positive parity excitation of the nucleon, the A^*(1440) Roper 
resonance, compared with the lowest lying odd parity excitation. In valence 



quark models with harmonic oscillator potentials, the J = ^ state naturally 

occurs below the N — 2, ^ excitation [6]. Without fine tuning of parameters, 
valence quark models tend to leave the Roper mass too high. Similar difficul- 
ties in the level orderings appear for the J^ = ^~^ Z\*(1600) and i"*" Z'*(1690) 
resonances, which have led to speculations that the Roper resonances may be 
more appropriately viewed as hybrid baryon states with explicitly excited glue 
field configurations [7] , rather than as radial excitations of a three-quark config- 
uration. Other scenarios portray the Roper resonances as "breathing modes" of 
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the ground states [8] , or states which can be described in terms of mcson-baryon 
dynamics alone [9]. 

Another challenge for spectroscopy is presented by the anomalously small 
mass of the odd parity /1(1405) hyperon. This is alternatively interpreted as a 
true three-quark state, or as a hadronic molecule arising from strong coupled 
channel effects between Stt, KN, ... states [10]. In fact, the role played by 
Goldstone bosons in baryon spectroscopy has received considerable attention 
recently [11-13]. 

Mass splittings between states within SU(3) quark- model multiplets provide 
another important motivation for studying the baryon spectrum. The dynamical 
origin of hyperfine splittings in quark models has traditionally been attributed 
to the colour- magnetic one gluon exchange mechanism [14]. On the other hand, 
there have been attempts recently to explain the hyperfine splittings and level 
orderings in terms of a spin-flavour interaction associated with the exchange of a 
pseudoscalar nonet of Goldstone bosons between quarks [15]. Understanding the 
mass splitting between the J^ = j iV*(1535) and | A^*(1520), or between 

the i Zi*(1620) and | Z\*(1700), can help identify the important mechanisms 
associated with the hyperfine interactions, and shed light on the spin-orbit force, 
which has been a central mystery in spectroscopy [16]. In valence quark models, 
the degeneracy between the N^ and A^| can be broken by a tensor force 
associated with mixing between the N"^ and N'^ representations of SU(3) [5], 
although this generally leaves the 7V| at a higher energy than the N^ .In 

contrast, a spin-orbit force is necessary to split the Z\| and Z\^ states. In the 
Goldstone boson exchange model [15], both of these pairs of states are degener- 
ate. 

Similarly, neither spin-flavour nor colour-magnetic interactions is able to ac- 
count for the mass sphtting between the yl(1405) and the J^ = j yl*(1520). A 
splitting between these can arise in constituent quark models with a spin-orbit 
interaction, but the required strength of the interaction leads to spurious mass 
splittings elsewhere [5,17]. In fact, model- independent analyses of baryons in 
the large Nc limit [18] have found that these mass splittings receive important 
contributions from operators that do not have a simple quark model interpreta- 
tion [18], such as those simultaneously coupling spin, isospin and orbital angular 
momentum, as anticipated by early QCD sum-rule analyses [16]. Of course, the 
coefficients of the various operators in such an analysis must be determined 
phenomenologically and guidance from lattice QCD is essential. 

The large number of "missing" resonances, predicted by the constituent quark 
model and its generalisations but not observed experimentally, presents another 
problem for spectroscopy. If these states do not exist, this may suggest that 
perhaps a quark-diquark picture (with fewer degrees of freedom) could afford 
a more efficient description. Alternatively, the missing states could simply have 
weak couplings to the TriV system [5] . 

Apart from these long-standing puzzles, interest in baryon spectroscopy has 
been further stimulated with the recent discovery of the exotic strangeness S = 
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+ 1 O^ pcntaquark [19], which has a minimuin 5 valence quark [uudds) content. 
Postulated long ago, the 0'^ appears to have eluded searches until now, and 
understanding its properties and internal structure has become a major challenge 
for spectroscopy. 

A consistent description of the baryon spectrum within low-energy models 
requires a satisfactory resolution of the old questions of spectroscopy, as well 
as insight into the new puzzles. On both of these fronts real progress may only 
come with the help of first-principle calculations of the spectrum in lattice QCD, 
which at present is the only method capable of determining hadron properties 
directly from the fundamental quark and gluon theory. Considerable progress 
has been made in the past few years in calculating hadronic properties in both 
quenched and full QCD, with the ground state masses now understood at the 
few percent level. 

Compared with simulations of hadron ground state properties, however, the 
calculation of the excited nucleon spectrum places particularly heavy demands 
on lattice spectroscopy. The excited nucleon states are expected to be large, 
with the size of a state expected to double with each increase in orbital angular 
momentum. Lattice studies of the excited nucleon spectrum therefore require 
large lattice volumes, with correspondingly large computational resources. Fur- 
thermore, the states are relatively massive, requiring a fine lattice spacing, at 
least in the temporal direction. Recent advances in computational capabilities 
and more efficient algorithms have enabled the first dedicated lattice QCD sim- 
ulations of the excited states of the nucleon to be undertaken. 

Of course, the calculation of the hadronic spectrum faces a formidable chal- 
lenge in describing excited state decays, which is relevant in both full and 
quenched QCD. Lattice studies of excited hadrons are possible because at the 
current unphysically large quark masses and finite volumes used in the simula- 
tions many excited states are stable [20]. Contact with experiment can be made 
via extrapolations incorporating the nonanalytic behaviour predicted by chiral 
effective field theory. 

The rest of this review is laid out as follows. After briefly reviewing in Sec- 
tion 2 the history of lattice calculations of the excited baryon spectrum, in 
Section 3 we outline the basic lattice techniques for extracting masses from cor- 
relation functions, including variational methods and Bayesian statistics. Sec- 
tion 4 describes the construction of a suitable basis of interpolating fields for 
any baryon we may wish to investigate. Section 5 outlines the simplest interpo- 
lating fields for spin-1/2 and spin-3/2 baryons, and how these can be used in the 
construction of lattice correlation functions. In Section 6, we present a survey of 
recent N* results for both positive and negative parity baryons. Finally, in Sec. 7 
we discuss conclusions of the existing studies, and speculate on future directions 
for the study of baryon spectra. Some technical aspects of the correlation matrix 
formalism for calculating masses, coupling strengths and optimal interpolating 
fields are described in the Appendix. 
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2 History of Lattice N* Calculations 

The history of excited baryons on the lattice is rather short, although recently 
there has been growing interest in finding new techniques to isolate excited 
baryons, motivated partly by the experimental N* programme at Jefferson Lab. 
An attempt to determine the first positive-parity excitation of the nucleon using 
variational methods was made by the UKQCD collaboration [21] . A more detailed 
analysis of the positive parity excitation of the nucleon was performed by Lein- 
weber [22] using Wilson fermions and an operator product expansion spectral 
ansatz. DeGrand and Hecht [23] used a wave function ansatz to access P-wave 
baryons, with Wilson fermions and relatively heavy quarks. Subsequently, Lee 
and Leinweber [24] introduced a parity projection technique to study the neg- 
ative parity ^ states using an 0{a?) tree- level tadpole- improved 0^34 quark 
action, and an O(a^) tree-level tadpole-improved gauge action. Following this, 
Lee [25] reported results using a D234 quark action with an improved gauge ac- 
tion on an anisotropic lattice to study the ^ and ^ excitations of the nucleon. 
An anisotropic lattice with an 0{a) improved quark action was also used by 
Nakajima et al. [26] to study excited states of octet and decuplet baryons. 

The RIKEN-BNL group [27] has stressed the importance of maintaining chi- 
ral symmetry on the lattice. At finite lattice spacing the Wilson fermion action 
is known to explicitly break the chiral symmetry of continuum QCD. A solu- 
tion to this problem is provided through the introduction of a fifth dimension, 
which allows chiral symmetry to be maintained even at non-zero lattice spacing 
[28]. The resulting domain wall fermion action was used in Ref. [27] to compute 
the masses of the N^ and N^ excited states. The analysis has recently been 

extended by studying the finite volume effects of the first N^ excited state 
extracted using maximum entropy methods [29] . 

A nonperturbatively O(a)-improved Sheikholeslami-Wohlert (SW) [30], or 
clover, fermion action was used by Richards et al. [31] to study the N^ and 
Z\| states. By appropriately choosing the coefficient of the improvement term, 
all 0{a) discretisation uncertainties can be removed, ensuring that continuum- 
like results are obtained at a finite lattice spacing. For more details on 0{a)- 
improvement, see the contribution by Zanotti et al. in this lecture series. 

The BGR collaboration [32] has been investigating the masses of the N^ 

and N^ excited states calculated with chirally improved (CI) and Fixed Point 
(FP) fermions. Both of these actions offer the advantage of improved chiral prop- 
erties over traditional Wilson-style fermions, without the cost associated with 
Ginsparg- Wilson fermions, which possess an exact analogue of chiral symmetry. 
The CSSM Lattice Collaboration has used an 0{a'^) improved gluon action 
and the C'(a)-improved Fat Link Irrelevant Clover (FLIC) fermion action [33] to 
perform a comprehensive study of the spectrum of positive and negative parity 
baryons [34,35]. Excited state masses in both the octet and decuplet multiplets 
have been computed, including the N^ , N^ , S^ , A^ , S^ , Z\i and 
A^ baryons. The formulation of FLIC fermions and a review of its scaling and 



Baryon Spectroscopy in Lattice QCD 5 

light quark mass phenomenology is reviewed by Zanotti et al. in this lecture 
series. 

Constrained-fitting methods originating from Bayesian priors have also re- 
cently been used by Lee et al. [36] to study the two lowest octet and decuplet 
positive and negative parity baryons, using overlap fermions with pion masses 
down to ^ 180 MeV. These authors have addressed the difficult issue of con- 
tamination of the first excited state of the nucleon with quenched rj'N artifacts 
[37]. While these authors claim to have observed the Roper in quenched QCD, it 
remains to be demonstrated that this conclusion is independent of the Bayesian- 
prior assumed in their analysis [22]. It would be interesting to examine these 
correlation functions using correlation matrix techniques or alternative Bayesian 
techniques such as the Maximum Entropy Method. 

As mentioned in Sec. 1, we are currently seeing intense interest in exotic 
pentaquark spectroscopy, and recently the first lattice studies of 5-quark states 
have appeared [38-40]. Here the unique advantage of lattice QCD can come 
to the fore, in predicting, for instance, the parity and spin of the lowest lying 
pentaquark state, which are as yet undetermined experimentally. The first two 
of these early studies have used the standard Wilson gauge and fermion action, 
on lattice sizes oi L = 1.2-2.7 fm [38], and L = 2.2 fm [39], while the third used 
domain wall fermions with a Wilson gauge action on a. L — 1.8 fm lattice. 

The main challenge here has been the construction of lattice operators for 
states with 4 quarks and an antiquark, which can be variously constructed in 
terms oi {qqq){qs) "nucleon-meson" [38] or {qq){qq)s "diquark-diquark-s" [39,40] 
operators. While two of these studies [38,39] appear to favor negative parity for 
the lowest lying pentaquark state, it is not clear that this state is a new resonant 
combination oi K N in S wave. The most recent study [40] reports an even 
parity ground state, and moreover identifies the Roper resonance as a pentaquark 
state. The relatively small volumes, limited set of interpolating fields, and naive 
linear extrapolations means that these results must be regarded as exploratory 
at present. New studies using the FLIC fermion action with several different 
interpolating fields are currently in progress [41]. 

We are clearly witnessing an exciting period in which the field of baryon 
spectroscopy on the lattice is beginning to flourish. 



3 Lattice Techniques 

3.1 Spectroscopy Recipe 

The computation of the spectrum of states in lattice QCD is in principle straight- 
forward. The building blocks are the quark propagators 

S%ix,y)^{0\i;:,ix)^^piy)\0), (1) 

computed on an ensemble of gauge configurations, composed of gauge field links 
Ufj_{x) ~ exp{mgAp(x)}. The calculation then proceeds as follows: 
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1. Choose an interpolating operator O that has a good overlap with P, the 
state of interest, 

(0 I O I P) ^ 0, 

and ideally a small overlap with other states having the same quantum num- 
bers. 

2. Form the time-sliced correlation function 

G(i,p) = ^(O(a:,t)Ot(0,0))e-'P-, 

which can be expressed by a Wick expansion in terms of the elemental quark 
propagators of Eq. (1). 

3. Insert a complete set of states between O and O^. The time-sliced sum puts 
the intermediate states at definite momentum, and one finds 



X P 



E 



{^^\0\P)^^_,EAv)t 



P ^Epip) 

where the sum over P includes the contributions from two-particle and higher 
states. 
4. Continue to Euclidean space t -^ —it, giving 

At large Euclidean times, the lightest state dominates the spectral sum in 
Eq. (2), and we can extract the ground state mass. The determination of this 
ground state mass for the states of lowest spin has been the benchmark calcula- 
tion of lattice QCD since its inception. However, our goal is to build up a more 
complete description of the baryon spectrum, ultimately determining not only 
the masses of some of the higher spin particles, but also the masses of the radial 
excitations. In the remainder of this section, we will address two issues that are 
crucial to attaining this goal: the application of variational techniques to isolate 
the higher excitations in Eq. (2), and the use of the appropriate statistical fitting 
techniques to reliably extract the energies of those excitations. In the following 
section, we will describe the construction of the nucleon interpolating operators. 

3.2 Fitting Techniques 

Variational Methods To confidently extract other than the lowest-lying state 
in Eq. (2), it is crucial to have more than a single interpolating operator Oi 
in order to appeal to variational methods to determine the spectrum of states 
[42,43]. Our aim is to measure every element of the matrix of correlators 

G,,{t)^{J20.{x,t)0]m. (3) 
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We now consider the eigenvalue equation 

G{t)u = X{t,to)G{to)u, (4) 

for some eigenfunction u of G~^{to)G{t), with tg fixed. For the sake of illustra- 
tion, we will consider a simple system with only two independent states, and a 
2x2 matrix of correlators. Then the eigenvalues of Eq. (4) satisfy 

A+(t,to)=e-(*-*''^^'' , 

A_(t,io)=e-(*-*''^^^ (5) 

and we have an exact separation between the energies of the two states, with 
coefficients growing exponentially with Iq. For the physical case of more than two 
states, there are exponential corrections arising from the states of higher energy. 
Ideally, to suppress the contribution of these higher energy states we would wish 
to choose to to have as large a value as possible. However, increasing to comes 
at the price of increasing statistical noise, and therefore we are generally obliged 
to take to close to the origin. 

In practice, it is usual to adopt a variation of the above method, and not to 
attempt to diagonalise the transfer matrix at each time slice, but rather to choose 
a matrix of eigenvectors V{to) for some to close to the source that diagonalises 
G{to)~^G{to + 1) . For the case of a 2 x 2 matrix with only two states, the 
diagonal elements of the matrix 

V{to + l)-^G{to)-'G{t)V{to + l) 

are indeed equivalent to the eigenvalues of G(io)^^G(io + 1). A more pedagogical 
discussion of these concepts can be found in the Appendix. 

The efficacy of this method relies on having a basis of correlators that can de- 
lineate the structure of the first few states, together with sufficiently high-quality 
statistics that the elements of the correlator matrix can be well determined. In 
the case of the glueball spectrum, the computational cost is dominated by the 
cost of generating the gauge configurations, and the overhead of measuring ex- 
tra correlators is negligible. Furthermore, we are able to improve the statistical 
quality of the data by using translational invariance to average over the position 
of the source coordinate in Eq. (2) . Variational techniques have been essential, 
and very successful, in the extraction of the glueball spectrum [44,45]. 

In the case of operators containing quark fields, there is generally a consider- 
able computational cost associated with the measurement of additional correla- 
tors. We will discuss the construction of such operators in the next section, but 
it is straightforward to see how this cost arises. Recall that the basic building 
blocks of spectroscopy are the quark propagators of Eq. (1), which satisfy 

M'J^ix,y)Sj,';iy,z)^S^H^^d,.. (6) 

where M is the fermion matrix of the fermion action tp Mip. The propagator is 
usually obtained by standard sparse matrix inversion methods, such as conjugate 
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gradient. These methods are apphcable to a fixed source vector on the right 
hand side of Eq. (6), and additional correlators require the inversion of M for 
additional source vectors. A corollary of this is that the matrix of correlators 
is obtained for fixed source coordinate, and therefore we cannot in general use 
translational invariance to average over the source coordinate in the manner of 
glueball calculations. 

Bayesian Statistics Given the computational cost of measuring additional 
correlators, it is important to extract the greatest possible information from the 
those that are measured. Furthermore, it is vital that the fitting procedure be as 
reliable as possible, in particular by ensuring that the extracted masses are not 
adversely affected by contamination from higher excitations. Historically, lattice 
calculations employed maximum likelyhood fits to correlators such as Eq. (2), 
using a single or at most two states in the spectral sum. An acceptable x^/d.o.f. 
would only be obtained if the fitted data were restricted to large temporal sepa- 
rations, in which one or two states did indeed dominate the spectral sum. Thus 
the fits ignored the data at small temporal separations, which we have already 
seen has the largest signal-to-noise ratio. 

It is natural to ask whether one can extract useful information from the data 
at small values of i, and a means of so doing is provided through the use of 
Bayesian Statistics. We will return to a further discussion of the the Bayesian 
approach in our survey of results. 

4 Interpolating Fields 

The variational techniques and fitting methods described above require a suitable 
basis of interpolating fields from which operators can be constructed which mimic 
the structure of each of the states to be extracted. In order to do so, we firstly 
have to consider the possible quantum numbers by which states are classified, and 
the extent to which they remain good quantum numbers in lattice calculations. 
We then have to consider the spatial structure of the states, and in particular 
their spatial extents, and whether we can construct interpolating operators that 
reflect the structure, through the use of, say, smeared interpolating fields. We 
will begin this section by detailing the quantum numbers with which states are 
classified, namely flavour and parity, and the degree to which they are good 
quantum numbers on the lattice. We will then proceed to discuss the use of 
spatially extended interpolating operators, known as smearing. Finally, we will 
describe the operators that have been employed thus far in lattice calculations, 
and discuss these in terms of continuum quark model operators. 

4.1 Continuum and Lattice Symmetries 

Baryon states are classified by their fiavour structure, either according to SU{2) 
(isospin), or SU{i) (strangeness), parity, and total spin. In nearly all lattice 
calculations of the hadron spectrum, exact isospin symmetry is imposed so the 
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TO„ — nid, and electromagnetic effects are ignored. Thus the flavour structure of 
baryon states composed of light (u, d, s) quarks is specified according to Total 
Isospin I, I3, and the strangeness S, and the naming of baryon states follows 
from this labeling, as detailed in Table 1 . 



Table 1. The table shows the flavour classification of baryon states constructed from 
light (u, d, s) quarks, together with a representative flavour structure in terms of three 
valence quarks for the case J3 = 7. 



/ 


s 


Baryon 


Flavour structure (is = /) 


1/2 





N 


^(1 udu)^ 1 duu)) 


3/2 





A 


uuu) 





-1 


A 


A={\ uds)- 1 dus)) 


1 


-1 


E 


1 uus) 


1/2 


-2 


S 


1 uss) 





-3 


n 


1 sss) 



The flavour structure is straightforward to implement in the calculation of 
the correlation functions by only including the appropriate contractions in the 
Wick expansion of Eq. (2). 

4.2 Angular Momentum and Lattice QCD 

More delicate is the identification of the spin of particles in a lattice calculation. 
For the discussion of the spectrum, we are principally interested in the study of 
states at rest, for which the relevant continuum symmetry for classifying states is 
their properties under rotations, described by the group SU{2). The irreducible 
representations of SU{2) are labeled by the total spin J, and the projection of 
spin along a particular axis, say J3, in the manner of isospin. We have already 
seen that we can impose exact isospin symmetry in our calculations, but the 
use of a hypercubic lattice has the consequence that we no longer have exact 
three-dimensional rotational symmetry, but rather the symmetries of the three- 
dimensional cubic group of the three-dimensional spatial lattice, the octahedral 
group O [46]. 

In contrast to the continuum rotation group, O is a finite group. It contains a 
total of 48 elements, and has a total of five single- valued irreducible representa- 
tions (IR), corresponding to states of integer spin, and only three double- valued 
IRs, corresponding to states of half-integer spin. Thus each lattice IR is an ad- 
mixture of different values of J, listed for the double-valued representations in 
Table 2 [47]. 

Thus in any lattice calculation at a fixed value of the lattice spacing, we 
will aim to construct operators transforming irreducibly according to one of the 
IRs of Table 2, and extract the spectrum of states within each of these IRs. We 
will only be able to identify the angular momentum of the various states when 
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Table 2. Irreducible representations of the octahedral group, listing the number of 
elements and the spin content. 



Representation 



Gi 
Ga 
H 



Dimension 



Spin (J) 



1/2,7/2,9/2,. 
5/2, 7/2, 11/2, 
3/2, 5/2, 7/2, . 



we look for commonality between the masses extracted from the various IRs in 
the approach to the continuum limit. Thus, for example, the continuum state of 
spin 5/2 has two degrees of freedom associated with G2^ and the remaining four 
degrees of freedom associated with H. 

So far, most lattice calculations have employed local, S'-wave propagators for 
the quarks, and as we shall see later, have employed operators transforming ac- 
cording to the IRs Gi, for spin ^ or H , for spin |, and have implicitly assumed 
that these states with these spins are indeed the lightest states in their respective 
representations, as observed in the physical hadron spectrum. The technology 
required to construct general baryon interpolating operators transforming irre- 
ducibly under O has now be developed [48,49], and a preliminary attempt at the 
spectrum using P-wave quark propagators has been made in ref. [50]. 

4.3 Parity 

The remaining symmetry that we must consider is parity, corresponding to the 
spatial-inversion operator /„ . Clearly, this is a good symmetry on our hypercubic 
lattice, yielding the point group Oh, with 96 elements, and an additional subscript 
g or u on our irreducible representations corresponding to positive- and negative- 
parity representations respectively. One subtlety arises when one considers the 
determination of the spectrum on a lattice with either periodic or anti-periodic 
boundary conditions in the temporal direction. In that case, for identical source 
and sink we have 

G{t) ^ Y.^o+{x,t)o+m 

X 

P+ P- 

where 0+ is an interpolating operator designed to isolate states of positive parity 
propagating in the forward time direction, Nt is the temporal extent of the 
lattice, Ap is the amplitude of a state P, and the subscripts -f and — denote 
contributions of states of positive or negative parity respectively. As discussed 
further in Sec. 5.1, a similar expression may be obtained for the opposite parity 
states. Thus we see that a complete delineation of the states on a periodic lattice 
only occurs as Nt -^ oo. The superposition of results from periodic and anti- 
periodic boundary conditions can be used to eliminate the second term of Eq. (7). 
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Alternatively, a fixed boundary condition can be used to eliminate the second 
term of Eq. (7) by preventing states from crossing the temporal boundary of the 
lattice. 



4.4 Smearing and Extended Interpolating Fields 

It is common to perform some smearing in the spatial dimensions at the source to 
increase the overlap of the interpolating operators with the ground states. Here 
we describe one such technique to do this: gauge-invariant Gaussian smearing 

[51]. 

The source-smearing technique [51] starts with a point source, 

^oli^^t) -- S'^'S^j^^^^okto (8) 

for source colour c, Dirac 7, position Xq and time to, and proceeds via the 
iterative scheme, 

x' 

where 

Repeating the procedure TV times gives the resulting fermion source, 

V'jv(a;,t)=^^^(a;,a;')V'o(a^',t). (9) 

x' 

The parameters N and a govern the size and shape of the smearing function. 
The propagator S is obtained from the smeared source by solving 

M^^^^^-^^C (10) 

for each colour, Dirac source c, 7 respectively of Eq. (8) via a standard matrix 
inverter such as the BiStabiHsed Conjugate Gradient algorithm [52]. 

5 Operators for Spin-| and Spin-| Baryons 

We will now see how the considerations described above apply to the construction 
of operators for the spin-i and spin-| baryons. Initially we will construct the 
operators based on analogy with the construction of the interpolating operators 
of the continuum, first for spin 2 and then for spin |. We will then briefly 
describe the construction of the operators that transform irreducibly under the 
lattice symmetries directly, without reference to the continuum discussion. 
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5.1 Spin-i Baryons: General Considerations 

Following standard notation, we define a two-point correlation function for a 
spin-i baryon B as 

GB{t,p) ^ Y. ^'^"'^ (^ Ixb(x)xb(0)| n) (11) 

where xb is a baryon interpolating field and where we have suppressed Dirac in- 
dices. All formalism for correlation functions and interpolating fields presented in 
this article is carried out using the Dirac representation of the Dirac 7-matrices. 
The choice of interpolating field xb is discussed in Section 5.2 below. The over- 
lap of the interpolating field xb with positive or negative parity states |-B ) 
is parameterised by a coupling strength \b± which is complex in general and 
which is defined by 



(^|XB(0)|B+,p,s)-AB+W^ws+fe5), (12) 



E 



B + 



{n\XB{0)\B-,p,s)=\B-J^75UB-{P,s), (13) 



where Mb± is the mass of the state B"^, Eb± = JM^^ + p^ is its energy, and 



Ub±{p,s) is a Dirac spinor with normalisation Ug±{p, s)u^±{p, s) — S"'^ . For 
large Euclidean time, the correlation function can be written as a sum of the 
lowest energy positive and negative parity contributions 

(^B(t,p)^ Ab+ TT^ e B+ + \g_ — e B , (14j 

when a fixed boundary condition in the time direction is used to remove back- 
ward propagating states. The positive and negative parity states are isolated by 
taking the trace of Gb with the operator F^ and F- respectively, where 

1 L . Mbt 



For p ~ 0, the energy Ebt — Mbt, so that Fi- — Up and the F:^ are parity 
projectors In this case, positive parity states propagate in the (1, 1) and (2, 2) 
elements of the Dirac 7-matrix of Eq. (14), while negative parity states propagate 
in the (3, 3) and (4, 4) elements, and the masses of B^ may be isolated. 

5.2 Spin-i Baryons: Interpolating Fields 

There are two types of interpolating fields which have commonly been used in 
the literature. The notation adopted here is similar to that of Rcf. [53]. To access 
the positive parity proton we use as interpolating fields 

Xi^(2;) = ^abc (ulix) C-fs db{x)) Uc{x) , (16) 
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and 



X2 i^) = ^atc {Ua {x) C db{x)) 75 Uc{x) 



(17) 



where the fields u, d are evaluated at Euclidean space-time point a;, C is the 
charge conjugation matrix, a, b and c are colour labels, and the superscript T 
denotes the transpose. These interpolating fields transform as spinors under a 
parity transformation. That is, if the quark fields qa{x) (q ^ u,d, ■ ■ •) transform 
as 

Vqa{x)V^ = +"/oqa{i) , (18) 

where x = (xq, —x), then 

Px^+(x)7't = +7oX^+(i). 
For convenience, we introduce the shorthand notation 

g{Sf„Sf,,Sf,) ^ e"^^6"''''^'|^;f'(x,0)tr Is'}'^' ^{x,0)S%' {x,0) 



^bb'T 



+ S^f^ix,0)S}yix,0)S%ix,0)\, 



(19) 



where S*?" (x, 0) are the quark propagators in the background link-field configu- 
ration U corresponding to flavours fi-3. This allows us to express the correlation 
functions in a compact form. The associated correlation function for Xi can be 
written as 



Gf+(t,p;r) = ( ^e-^P-^tr l-F G (Su, CSdC-\ Su 



(20) 

where (• • •) is the ensemble average over the link fields, F is the /± projection 
operator from Eq. (15), and C — C75. For ease of notation, we will drop the 
angled brackets, (• • •), and all the following correlation functions will be under- 
stood to be ensemble averages. For the X2 interpolating field, one can similarly 
write 

Glt{t,p; n^Y. ^''""'^^ [-^ ^ (75^«75, CSdC'\ 755«75)] , (21) 

X 

while the interference terms from these two interpolating fields are given by 

G\+{t,p-r) = ^e-^P-tr -r Ig (5^75, CSdC-\ Sul^)\ , (22) 

Gl+{t,p-r) = ^e-'P-tr -r Ig (755^, CSdC-\ 75^«)| . (23) 

The neutron interpolating field is obtained via the exchange u ^^ d, and the 
strangeness -2, S interpolating field by replacing the doubly represented u or d 
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quark fields in Eqs. (16) and (17) by s quark fields. The S and S interpolators 
are discussed in detail below. 

As pointed out in Ref . [22] , because of the Dirac structure of the "diquark" 
in the parentheses in Eq. (16), in the Dirac representation the field Xi involves 
both products of upper x upper x upper and lower x lower x upper components 
of spinors for positive parity baryons, so that in the nonrelativistic limit x^ = 
C(l). Here upper and lower refer to the large and small spinor components 
in the standard Dirac representation of the 7-matrices. Furthermore, since the 
"diquark" couples to total spin 0, one expects an attractive force between the 
two quarks, and hence better overlap with a lower energy state than with a state 
in which two quarks do not couple to spin 0. 

The X2 interpolating field, on the other hand, is known to have little 
overlap with the nucleon ground state [22,54]. Inspection of the structure of 
the Dirac 7-matrices in Eq. (17) reveals that it involves only products of up- 
per X lower X lower components for positive parity baryons, so that xf = 
OijP' / E^) vanishes in the nonrelativistic limit. As a result of the mixing of up- 
per and lower components, the "diquark" term contains a factor cr • p, meaning 
that the quarks no longer couple to spin 0, but are in a relative i = 1 state. 
One expects therefore that two-point correlation functions constructed from the 
interpolating field X2 ^^'^ dominated by larger mass states than those arising 
from Xi at early Euclidean times. 

While the masses of negative parity baryons are obtained directly from the 
(positive parity) interpolating fields in Eqs. (16) and (17) by using the parity 
projectors /^±, it is instructive nevertheless to examine the general properties of 
the negative parity interpolating fields. Interpolating fields with strong overlap 
with the negative parity proton can be constructed by multiplying the previous 
positive parity interpolating fields by 75, Xi^ — 75 Xi2- ^^ contrast to the 
positive parity case, both the interpolating fields Xi~ a-nd X2^ mix upper and 
lower components, and consequently both Xi~ and X2 are 0{p/E). 

Physically, two nearby J^ = ^ states are observed in the excited nucleon 
spectrum. In simple quark models, the splitting of these two orthogonal states is 
largely attributed to the extent to which the wave function is composed of scalar 
diquark configurations. It is reasonable to expect Xi~ to have better overlap with 
scalar diquark dominated states, and thus provide a lower effective mass in the 
moderately large Euclidean time regime explored in lattice simulations. If the 
effective mass associated with the X2~ correlator is larger, then this would be 
evidence of significant overlap of X2~ with the higher lying N^~ states. In this 
event, a correlation matrix analysis (see Appendix) can be used to isolate these 
two states [32,34]. 

Interpolating fields for the other members of the flavour SU(3) octet are 
constructed along similar lines. For the positive parity 17° hyperon one uses [53] 

Xf (a;) = ^eabc{ (ulix) C75 Sb{x)) dc{x) + (dlix) C75 Sfc(x)) u,(x)| , (24) 



Baryon Spectroscopy in Lattice QCD 15 

xf{x)^-^<^abc[{ua{x)Csbix)) -f5 dc{x) + (d^ (x) C Sb{x)) -f5Ucix)j. (25) 

Interpolating fields for accessing other charge states of S are obtained by d ^ u 
OTU ^ d, producing correlation functions analogous to those in Eqs. (20) through 
(23). Note that Xi transforms as a triplet under SU(2) isospin. An SU(2) singlet 
interpolating field can be constructed by replacing " + " — > " — " in Eqs. (24) 
and (25). For the SU(3) octet A interpolating field (denoted by "/l^"), one has 

Xi (x) = ^^abc {2 (ua (x) C75 db{x)) Sc{x) + {ul{x) C75 Sb{x)) dc{x) 

- (dlix) C75 Sb{x)) u,{x)} , (26) 
xf{x) = — =eQ6c|2(wf(a;)C4(x)) 75 Sc(x) + (^^(x) Csb(x)) -fbddx) 

- {dl{x)Csb{x)) j5Uc{x)}, (27) 
which leads to the correlation function 

G^;(i,p;r) = i^e-P- 

as 

-1(20 (^Ss, Sd, cSuC-^)+2g(^Ss, Su, cSdC-') 

+2Q iSd, Ss, CSuC \ +2Q [Su, &, CSdC 

- g (Sd. Su, cSsC-') - g (5„, Sd, cSsC-^J | , (28) 

and similarly for the correlation functions G22 , G^g and Gifi ■ 

The interpolating fields for the SU(3) flavour singlet (denoted by "A^") are 
given by [53] 

Xi (x) = -2 eabc {- {ua{x) G75 db{x)) Sc{x) + (u^(x) G75 Sb{x)) dc{x) 

- (dlix) Cj, Sb{x)) u,{x)} , (29) 

X2 (x) = -2 eabc\^ - {ul (x) C db{x)) 75Sc(x) + {u^ (x) C Sbix)) 75 4(x) 

-(C(x)Gsb(x)) 75Wc(x)} , (30) 

where the last two terms are common to both Xi a-nd Xi ■ The correlation 
function resulting from this field involves quite a few terms, 

GtUt,p;r) = 



xtr 



J2 e^'P^'tT 



-r \ 755^' CSf' ^G-i Sf-ir, + -i^sf CSf ^G-i 5f' 75 
+ j5Sr' CS^/^C-^ Sfj5 + l5Sf CS^/^C-^ 5f' 75 
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- 75^r'75 tr 



- 15 ST 75 tr 
l5Sfj5 tr 



obb' /~i Qcc T /~i—l 

obb' /^ Qcc T /~i—l 



Qbb r~i Qcc 1 ri 



' Tn-1 



(31) 



In order to test the extent to which SU(3) flavour symmetry is vahd in the 
baryon spectrum, one can construct another A interpolating field composed of 
the terms common to A^ and A^, which does not make any assumptions about 
the SU(3) flavour symmetry properties of A. We define 

xf{x) = --^<^abc\^ {ul{x) C75 Sb{x)) dc{x) - {dl{x) C75 Sh{x)) Uc{x)^ , (32) 
X2'i^) = -7^<^abc\^ {ul{x)C Sb{x)) 75 4(a;) - {dl{x)C Sh{x)) 75 Mc(a;)| , (33) 

to be our "common" interpolating fields which are the isosinglet analogue of 
Xi and X2 i'^ Eqs. (24) and (25). Such interpolating fields may be useful in 
determining the nature of the /I* (1405) resonance, as they allow for mixing 
between singlet and octet states induced by SU(3) flavour symmetry breaking. To 
appreciate the structure of the "common" correlation function, one can introduce 
the function 



G(St\,Sf2,Sfy 



abc a b c 



bb' T ( 



Sf^ix,0)tT S'}^ ' {x,0)S% ix,0) 



S'}f{x,0)S'}'^'^{x,0)S%{x,0) 



(34) 



which is recognised as G in Eq. (19) with the relative sign of the two terms 
changed. With this notation, the correlation function corresponding to the Xi 
interpolating field is 

G^f (t,p;r) = i^e-^P-^r -r {g (^Sd, CSsC-\ Su 



G [Su, CSgC , Sdj > 



(35) 



and similarly for the correlation functions involving the X2 interpolating field. 



5.3 Spin-2 Baryons 

In this section we extend the discussion from the previous section to the spin-| 
baryon sector. The mass of a spin-| baryon on the lattice is obtained from the 
two-point correlation function Gp,y [35,55], 

G^4t,p;r)=tv[rG^At,P)], (36) 
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where 

5::Sit,p) = E ^^'"'^ (^1 T (x"(^) xf (0)) \f2) , (37) 

X 

and Xu is a spin-| interpolating field, _r is a matrix in Dirac space with a, /3 
Dirac indices, and /x, i^ Lorentz indices. 

The interpolating field operator used in the literature for accessing the isospin- 
i, spin-|, positive parity (charge +1) state is [35,56,57] 



^p ^ ,abc ^^Ta^^^ c^^^. ^6(^)) (^^^^ _ 1^^^^^ ^^^^c(^) 



(38) 



As pointed out in Section 5.1, all discussions of interpolating fields are carried out 
using the Dirac representation of the 7-matrices. This exact isospin-^ interpolat- 
ing field has overlap with both spin-| and spin-i states and with states of both 
parities. The resulting correlation function will thus require both spin and parity 
projection. The charge neutral interpolating field is obtained by interchanging 
u ^-f d. This interpolating field transforms as a Rarita-Schwinger operator un- 
der parity transformations. That is, if the quark field operators transform as in 
Eq. (18), then 

and similarly for the Rarita-Schwinger operator 

ru^{x)V^ = +7oWm(^) > (39) 

which will be discussed later. 

The computational cost of evaluating each of the Lorentz combinations in 
Eq. (38) is relatively high — about 100 times that for the ground state nucleon 
[24] . Consequently, in order to maximise statistics it is common to consider only 
the leading term proportional to g^i, in the interpolating field, 

Xl — 6""= (u^'^ix) C757M d'ix)) j,u%x) . (40) 

This is sufficient since we will in either case need to perform a spin-| projection. 
In order to show that the interpolating field defined in Eq. (40) has isospin i, 
we first consider the standard proton interpolating field given in Eq. (16) which 
we know to have isospin 2 • Applying the isospin raising operator, /"*", on x^, one 
finds 

j+^P ^ e«''=(u^"C'75u'')u" 

= 0. 

Similarly, for the interpolating field defined in Eq. (40), one has 
I+Xl - e"'^(^^^"C757M"')75^^^ 
= e"^^(«^"C757M"')'^75"^ 
= -e"''^(u^"C757M«')75W^ 
= 0, 
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where we have used the representation-independent identities C^i^C^^ — — 7J, 
C^^C^^ = 7I", and the identities which hold in the Dirac representation: C"^ = 
C''' = C^^ ~ — C with C = 17270 and 7J = 75. 

We note that Xu, corresponding to yf in Eq. (40), is 



so that 



xP=e«''-r75(d'7M75Cu^^) , (41) 



Xlxl = e'^^^6«'^'^'(u^'^[C757M]a^4)75<<'75(4 [lul,C]p.o..ul?') 



—> 755'«75tr 'f^Sulb {C"ff_iSd-luC) 

IbSulb (Cjf^Sdj^C) 75S'«75 , (42) 

where the last line is the result after performing the Grassman integration over 
the quark fields with the quark fields being replaced by all possible pairwise 
contractions. 

In deriving the A interpolating fields, it is simplest to begin with the state 
containing only valence u quarks, namely the Z\+"'". The interpolating field for 
the Z\++ resonance is given by [57] 

xf^ix) = £«"= {u^^ix) Cj, u\x)) u^{x) , (43) 

which transforms as a pseudovector under parity. The interpolating field for a 
Z\+ state can be similarly constructed [55], 

xf{x) = -^ e^'^^ [2 {u^'^{x)C-i,d\x)) u%x) + {u^^{x) Cj,u''{x)) d^x)] . 

(44) 
Interpolating fields for other decuplet baryons are obtained by appropriate sub- 
stitutions of M, d -^ u, d OT s fields. 

To project a pure spin-| state from the correlation function G^^, one needs 
to use an appropriate spin-^ projection operator [58], 

Pjli^ (P) = .9m>- - 3 7m7<. - ^{l-Pl^iPi^+P^ili^l-p) ■ (45) 

The corresponding spin-^ state can be projected by applying the projection 
operator 

P'^i' = 9^u - Pl'J" ■ (46) 

To use this operator and retain all Lorentz components, one must calculate the 
full 4x4 matrix in Dirac and Lorentz space. However, to extract a mass, only one 
pair of Lorentz indices is needed, reducing the amount of calculations required by 
a factor of four. The results from Ref. [35] which are summarised in Section 6.4 
are calculated from the third row of the Lorentz matrix and using the projection 

4 
G^3= E ^3^5^"^ ^.'3, (47) 
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to extract the desired spin states, s == i or |. Following spin projection, the 
resulting correlation function, G'^^, still contains positive and negative parity 
states. 

The interpolating field operators defined in Eqs. (38) and (40) have overlap 
with both spin- 1 and spin- 2 states with positive and negative parity. The field 
Xp, transforms as a pseudovector under parity, as does the Rarita-Schwinger 
spinor, u^. Thus the overlap of Xp, with baryons can be expressed as 



ri + (^ .\\- \_,., l¥^ 



2+ 



{f2\x^\m + ip,s))=\y2+J-^ u^{p,s) , (48) 

V ^3/2+ 



3 Mr,/2- 

{f2\x^\m-{p,s)) = A3/2-W-^^ l5U^ip,s) , (49) 

V ^3/2- 



M1/2+ 



{f2\x^,\N-^ + ip,s)) ^ (ai/2+PM + /3i/2+7M)l IF^ 75"(p,s) , (50) 



E 



1/2 



+ 



i_ Mi/2- 

{n\x^,\N^ {j),s)) = (Qi/2-Pp + /?i/2-7m)i/ f "(P' *) ' (^1) 

y -f^i/2- 

where the factors A^ , oig , (3b denote the coupling strengths of the interpolating 
field x^l, to the baryon B. For the expressions in Eqs. (50) and (51), we note that 
the spatial components of momentum, pi, transform as vectors under parity and 
commute with 70, whereas the 7^ do not change sign under parity but anti- 
commute with 79. Hence the right-hand-side of Eq. (50) also transforms as a 
pseudovector under parity in accord with X/^- 
Similar expressions can also be written for x^. 



{m+{p,s)\xAn) = A*/2+W^2Z2l u,{p,s) , (52) 

V -C/3/2+ 



3 M-il-)- 

{m-{jp,s)\x^Q) = -XI .-^ S^(p,s)75 , (53) 

' V ^3/2- 



(7V2 + (p,s)|x^|/?) = -\ ^r^ "(P,s)75(aJ/2+PM + /?i/2+7M) ^ (54) 
V ^1/2+ 



1 / JKIl /9- 

(m-{p,s)\x,.\n) - J^^ S(p,s)(«t/2-PM + /3i*/2-7M) • (55) 



Note that we are assuming identical sinks and sources in these equations. How- 
ever, calculations often use a smeared source and a point sink in which case 
A*, a* and (3* are no longer complex conjugates of A, a and /3 but are instead 
replaced by A, a and (3. 

We are now in a position to find the form of Eq. (37) after we insert a complete 
set of intermediate states {\B{p,s))}. The contribution to Eq. (37) from each 
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intermediate state considered is given by 

{Q\x^Ni + {p,s)){Ni+{p,s)\xA^) 

M3/24 



=+A3/2+A3/2+ — u^,{p,s)u^{j),s) 

-C^3/2+ 



_ , Y M3/2+ (7 • P + M^/2+ ) J 1 2pf,p^ pf,j^ - p^^jfj, 

— A3/2+'\3/2+ -^ :^TT S g^iu - -tTmT!/ - XT72 1 T77 

{n\x^\Ni-{p,s)){Nl-{p,s)\x,\0) 

- M3/2- 
=-h/2-h/2- -e; l5Uf,{p,s)u^{p,s)jr, 

1^3/2- 
, Y ^3/2- il-p- M3/2-) \ 1 2p^,p^ Pf,J^ - p^7^ 

— A3/2-A3/2- ^3/^_ 2M3/2- V"^"^ r^'"" uq^ 3M3/2- 

(r2|xp|ivH(p,s))(7vH(p,s)|x.|^) 

== ^^ (ai/2+Pf, + Pi/2+7m)75 77T7 l5[ai/2+Pi^ + Pl/2+liy) > 

-C/1/2+ ^^^1/2+ 

Ml/2- N 7 ■ P + -^1/2- ,_ — 

=T^ (ai/2-PM + Pl/2-lf.) 77J7 {ai/2-P^ + (il/2-lu) ■ 

^-/2+ ^-'"1/2- 

To reduce computational expense, we consider the specific case when 11 = v = 3 
and in order to extract masses we require p = (0, 0, 0). In this case we have the 
simple expressions 



{0\X3\N^^{p,s)){Ni+{p,s)\xs\0) = A3/2+A3/2+^ ^l^^^l^l^^ , (56) 

(^|X3|^^-(p,5))(^i-(p,5)|x3|^> = A3/2-A3/2-^ ( ^"^'2^3'-^^'^'' ) '^^^^ 



{n\X3\N^+ip,s)){Ni+ip,s)\x3\0) - -/3i/2+/3i/2+7375 ^"^^m ^ ^'^'^ ^^73 

ZMI/2+ 

~^R 75 7oA^i/2+ + M1/2+ , . 

- +P1/2+P1/2+ ^77 . (oSJ 

(^IX3|iV^-(p,.))(JVi-(p,.)|X3|^)=/3i/2-^i/2-73 ^°^^y;."^^^^'' 73 

= +/9i/2- /3i/2- ^TT • (59) 

Therefore, in an analogous procedure to that used in Section. 5.1, where a fixed 
boundary condition is used in the time direction, positive and negative parity 
states are obtained by taking the trace of the spin-projected correlation function, 
GI3, in Eq. (47) with the operator F = r±, 

G^f^tr[r±Gy, (60) 
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where in this case (cf. Eq. (15)) 

^±-^(l±74). (61) 

The positive parity states propagate in the (1,1) and (2,2) elements of the Dirac 
7-matrix, while negative parity states propagate in the (3,3) and (4,4) elements 
for both spin-2 and spin-| projected states. A similar treatment can be carried 
out for the A interpolating fields but is left as an exercise for the interested 
reader. 

5.4 Octahedral Group Irreducible Representations 

The procedure outlined above largely follows the continuum construction of in- 
terpolating operators used in, say, QCD sum- rule calculations. An alternative 
approach is to construct operators that lie in the irreducible representations 
(IRs) of the cubic group Oh of the lattice directly [48,49]. To illustrate the latter 
method, we consider the construction of the / = 1/2, I^ = 1/2 nuclcon interpo- 
lating operator, following the discussion of Ref. [49]. 

The starting point is the formation of a set of elemental baryon operators 
that are gauge invariant, and have the correct isospin, or flavour, properties. In 
the case of point-like quark fields, these are given by 

'Pc.p-r = e'^'^«d>^ - d»^) , (62) 

where the coordinate of the quark fields are suppressed, and a, /3, 7 are spinor 
indices. Examination of Eq. (62) reveals the constraints 'Pa/ij + ^7/3a = and 
^ai3j + ^i3ja + ^70/3 = SO that thcrc are only 20 independent operators. 

The central result needed to find operators transforming irreducibly under 
the IRs of Table 2 is the projection formula 

Bf'^{x)^^J2DiTimRBf{x)UJ, , (63) 

where A refers to an Oh IR, A is the IR row, go^ is the number of elements 
in Oh, dA is the dimension of the A IR, Dmn{R) is a /I representation ma- 
trix corresponding to group element R, and Ur is the quantum operator which 
implements the symmetry operations; the temporal argument is suppressed. Ap- 
plication of the formula requires explicit representation matrices for each group 
element. These are applied to the 20 linearly independent operators above, from 
which 20 linearly independent operators that transform irreducibly under Oh 
are identified. The procedure is Dirac-basis dependent, and linearly independent 
operators are shown in Table 3, using the DeGrand-Rossi representation of the 
Dirac 7- matrices [49]. 

There are three embeddings of Gi and a single embedding of H for both even 
parity {g) and odd parity (u); the two helicities of spin 1/2 correspond to the two 
rows of the two-dimensional representation Gig(iti) while the four helicities of 
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spin 3/2 correspond to the four rows of the four-dimensional representation H. 
The three embeddings of Gi correspond to some hnear combinations of the op- 
erators of Eqs. (16)-(17) and the spin-1/2 projection of Eq. (38); the embedding 
of H corresponds to the spin-3/2 projection. 

Whilst the results in Table 3 may be less intuitive than their counterparts 
derived earlier, the procedure can easily extended to more complicated operators, 
including those with excited glue, or to pentaquark operators. Furthermore, we 
have seen explicitly that there are three independent spin-1/2 operators. 



Table 3. Combinations of the operators #0/37 in Eq. (63) which transform irreducibly 
under Oh for the DeGrand- Rossi representation of the Dirac 7-matrices, employed by 
LHPC and MILC Collaborations (see Ref [49]). 



Irrep 


Row 






Operators 


Gig 

Gig 


1 
2 






$112 + $334 
-$221 - $443 


Gig 

Gig 


1 
2 






$123 - $213 + $314 
$124 - $214 + $324 


Gig 

Gig 


1 
2 


2<?^114 
— 2^223 


+ 2#332 - ^123 - ^213 + 2<?134 " ^314 
- 2<f44i + $124 - $214 - 2$234 + $324 


Giu 

Giu 


1 
2 






$112 - $334 
-$221 + $443 


Giu 
Giu 


1 
2 






$123 - $213 - $314 
$124 - $214 - $324 


Giu 

Giu 


1 
2 


2<?ill4 
-2<l>223 


- 2$332 - $123 - $213 - 2$134 + $314 
+ 2$441 + $124 + $214 + 2$234 - $324 



tig 1 \/3($ll3 -f $33lJ 

Hg 2 $114 + $332 + $123 + $213 - 2$134 + $314 

Hg 3 $223 - $441 + $124 + $214 - 2$234 "f $324 

Hg 4 \/3($224 + $442) 

ITu i V3($ll3-$33l) 

Hu 2 $114 - $332 -I- $123 + $213 + 2$134 - $314 

Hu 3 $223 - $441 -I- $124 -|- $214 + 2$234 - $324 

Hu 4 \/3($224 - $442) 



6 Survey of Results 

Now that we have provided a detailed description of the procedure for the cal- 
culation of excited baryons on the lattice, we will now present an overview of 
recent lattice calculations of the excited baryon spectrum. The emphasis of most 
calculations has been the demonstration that the excited nucleon spectrum is 
indeed accessible to lattice calculations, and most of the calculations share sev- 
eral features. Namely, they are in the quenched approximation to QCD, are 
obtained on lattice volumes of 2 — 3 fm, employ pseudoscalar masses of around 
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500 MeV, and make fairly limited investigations of the systematic uncertainties 
on the calculations. 

The most comprehensive study of the excited nucleon spectrum has been ob- 
tained using FLIC fermions, detailed in Refs. [34,35] , and we will describe these 
calculations at length in these lectures. For some of the lowest-lying resonances, 
such as the Roper resonance and the odd-parity partner of the nucleon, there 
have been studies of the systematic uncertainties in the calculations, in partic- 
ular arising from the finite lattice spacing, and the need to perform a chiral 
extrapolation [13]. 

For completeness, we will first briefly describe the gauge and fermion actions 
used in the FLIC fermion analysis. Additional details of the simulations can be 
found in Refs. [33-35]. 

6.1 Lattice Actions for FLIC Calculation 

For the gauge fields, a mean- field improved plaquette plus rectangle action is 
used. The simulations are performed on a 16"^ x 32 lattice at (3 = 4.60, which 
corresponds to a lattice spacing of a = 0.122(2) fm set by a string tension analysis 
incorporating the lattice Coulomb potential [59] with ^/a — 440 MeV. For the 
quark fields, a Fat-Link Irrelevant Clover (FLIC) [33] action is implemented. A 
detailed description of the FLIC fermion action can be found in the contribution 
by Zanotti et al. in this series. 

A fixed boundary condition in the time direction is used for the fermions 
by setting Ut{x, Nt) = V a; in the hopping terms of the fermion action, with 
periodic boundary conditions imposed in the spatial directions. Gauge-invariant 
Gaussian smearing [51] in the spatial dimensions is applied at the source to 
increase the overlap of the interpolating operators with the ground states. 

Five masses are used in the calculations [33] and the strange quark mass is 
taken to be the second heaviest quark mass in each case. The analysis is based on 
a sample of 400 configurations, and the error analysis is performed by a third- 
order, single-elimination jackknife, with the x^ per degree of freedom (A'df) 
obtained via covariance matrix fits. 

6.2 Spin-i Baryons 

In Fig. 1 we show the nucleon and N*{^ ) masses as a function of the pseu- 
doscalar meson mass squared, m^. The filled squares indicate the results for the 
FLIC action, and the stars for the Wilson action [34] (the Wilson points are 
obtained from a sample of 50 configurations). We note here that the spatial size 
of the lattice is L = 1.95 fm and that from the values of rn^ given in Table 1 of 
Ref. [34] we have tti^X > 5.52. 

For comparison, we also show results from earlier simulations with domain 
wall fermions (DWF) [28] (filled triangles), a nonperturbatively (NP) improved 
clover action on anisotropic lattices at three different lattice spacings [26] (open 
diamonds), an NP improved clover action at /? = 6.2 [31] (open squares, open 
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Fig. 1. Masses of the nucleon (A*') and the lowest J^ — i excitation ('W"). The 
FLIC and Wilson results are from Ref [34], the DWF [28], Fixed Point [32], Chirally 
Improved (CI) Wilson [32], NP improved clover [31] and NP improved anisotropic clover 
[26] results are shown for comparison. The empirical nucleon and low lying N*{^ ) 
masses are indicated by the asterisks along the ordinate. 



circles and filled diamonds), and results using Fixed Point (FP) (crosses) and 
Chirally Improved (CI) (open inverted triangles) [32] fermion actions. The scat- 
ter of the different NP improved results is due to different source smearing and 
volume effects: the open squares are obtained by using fuzzed sources and local 
sinks, the open circles use Jaeobi smearing at both the source and sink, while 
the filled diamonds, which extend to smaller quark masses, are obtained from 
a larger lattice (32^ x 64) using Jaeobi smearing. The empirical masses of the 
nucleon and the three lowest ^ excitations are indicated by the asterisks along 
the ordinate. In an unquenched calculation, the results may shift by the order 
of 10% compared with a quenched calculation [12]. 

There is excellent agreement between the different improved actions for the 
nucleon mass, in particular between the FLIC [34], DWF [28], NP improved 
clover [26,31], FP and CI [32] results. On the other hand, the Wilson results lie 
systematically low in comparison to these due to the large 0{a) errors in this 
action [33]. A similar pattern is repeated for the N*{^ ) masses. Namely, the 
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Fig. 2. The masses of the lowest-lying positive- and negative-parity baryons in units 
of Vq^ versus a^ in units of rg, where ro is Sommer's scale [60]. The lines are linear fits 
in a^ /ro to the positive- and negative-parity baryon masses; different plotting symbols 
correspond to different types of smearing. Also shown are the physical values. 



FLIC, DWF, NP improved clover, FP and CI masses are in good agreement with 
each other, while the Wilson results again lie systematically lower. A mass split- 
ting of around 400 MeV is clearly visible between the N and N* for all actions, 
including the Wilson action, despite its poor chiral properties. Furthermore, the 
trend of the N*{^ ) data with decreasing tti^t is consistent with the approach 
to the mass of the lowest-lying physical negative parity N* states. 

In the case of the NP-improved clover fermion action, with 0{a^) discreti- 
sation errors, the calculation has been performed at three values of the lattice 
spacing, enabling a continuum extrapolation to be performed [31]. This is shown 
in Fig. 2, although it should be noted that a simple linear chiral extrapolation 
was performed in this calculation. There is a suggestion of a somewhat larger lat- 
tice spacing dependence for the higher excited resonance, emphasising the need 
to perform a careful analysis of systematic uncertainties in future calculations. 

Figure 3 shows the mass of the J^ = i states (the excited state is denoted 
by 'W(l/2+)"). As is long known, the positive parity X2 interpolating field does 
not have good overlap with the nucleon ground state [22] and a correlation matrix 
analysis confirms this result [34], as discussed below. It has been speculated that 
X2 niay have overlap with the lowest ^ excited state, the iV*(1440) Roper 
resonance [28]. In addition to the FLIC and Wilson results from the present 
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analysis, we also show in Fig. 3 the DWF results [28], and those from an earlier 
analysis with Wilson fermions together with the operator product expansion 
(OPE) [22]. The physical values of the lowest three i excitations of the nucleon 
are indicated by the asterisks. 

The most striking feature of the data is the relatively large excitation en- 
ergy of the N'{^ ), some 1 GeV above the nucleon. There is little evidence, 
therefore, that this state is the N*{\AAQ) Roper resonance. While it is possible 
that the Roper resonance may have a strong nonlinear dependence on the quark 
mass at rn^ ^ 0.2 GeV^, arising from, for example, pion loop corrections, it is 
unlikely that this behaviour would be so dramatically different from that of the 
A^*(1535) so as to reverse the level ordering obtained from the lattice. A more 
likely explanation is that the X2 interpolating field does not have good overlap 
with either the nucleon or the iV*(1440), but rather (a combination of) excited 

I state(s). 

Recall that in a constituent quark model in a harmonic oscillator basis, the 
mass of the lowest mass state with the Roper quantum numbers is higher than 
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the lowest P-wave excitation. It seems that neither the lattice data (at large 
quark masses and with our interpolating fields) nor the constituent quark model 
have good overlap with the Roper resonance. Better overlap with the Roper is 
likely to require more exotic interpolating fields. 

As mentioned in Section 2, Lee et al. [36] have performed a calculation using 
overlap fermions with pion masses down to ~ 180 MeV. Using new constrained 
curve fitting techniques, they extract excited states from a single correlation 
function calculated with the standard nucleon interpolating field in Eq. (16). 
The results from this calculation exhibit a dramatic drop in the mass of the first 
excited i state of the nucleon at light pion masses, reversing the level ordering 

of the first ^ and ^ excited states. It is important, however, that this result 
be shown to be independent of the constrained curve fitting techniques adopted 
for this analysis. A correlation matrix analysis involving several operators would 
shed considerable light on this issue. 

Recently there has also been speculation that the Roper resonance suffers 
from large finite volume errors [29]. To study this issue, the authors of Ref. [29] 
calculate the nucleon and its first positive and negative parity excitations on 
three different lattice volumes {La = 1.5, 2.2 and 3.0 fm). Using Maximum 
Entropy Methods, they find that on large volume lattices {^ 3.0 fm) the mass 
of the i excited nucleon state is reduced. A similar analysis remains to be 
performed for the first ^ nucleon state obtained from the same fermion action. 
At present Wilson fermion scaling violations allow the ^ excited nucleon state 

to sit lower than the first ^ nucleon state obtained from the improved DWF 
action. It is essential to compare the masses of these states using the same 
fermion actions to remove systematic errors such as these. Similarly, a correlation 
matrix analysis involving several operators remains desirable. 

The BGR Collaboration has performed a calculation of the excited nucleon 
spectrum on two lattice volumes {La ~ 1.8 and 2.4 fm) using Fixed Point and 

Chirally Improved Wilson fermions [32]. Using three different operators to cre- 

1 + 1 ~ 

ate the nucleon states, the lowest two states in both the ^ and ^ channels 

are extracted using a 3 x 3 correlation matrix. The eigenvectors and eigenvalues 
determined from the correlation matrix analysis should correspond to the solu- 
tions of Eq. (4). It is unclear, however, whether the constraints which have been 
employed in Ref. [32] are physical, and in particular whether the eigenvectors 
correspond to the optimal projections onto physical states defined by the eigen- 
vectors u of Eq. (4). With this caveat, a splitting between the two ^ states 

is identified and the excited ^ state is observed to sit above the i states, in 
agreement with Ref. [34] (see next section). The results on two volumes do not 
exhibit any volume dependence. 

In the next section we present results from a 2 x 2 correlation matrix analysis 
using FLIC fermions on a single lattice volume. 
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Fig. 4. Masses of the J^ — | and i nucleon states, for the FLIC action [34]. The 
positive (negative) parity states are labeled A''i [N^) and N2 {N2). The results from 
the projection of the correlation matrix are shown by the filled symbols, whereas the 
results from the standard fits to the XiXi and X2X2 correlation functions are shown by 
the open symbols (offset to the right for clarity). Empirical masses of the low lying i 
states are indicated by the asterisks. 



6.3 Resolving the resonances 

The mass splitting between the two lightest N*{^~) states (iV*(1535) & iV*(1650)) 
can be studied by considering the odd parity content of the xi ^-nd X2 interpolat- 
ing fields in Eqs. (16) and (17). Recall that the "diquarks" in xi and X2 couple 
differently to spin, so that even though the correlation functions built up from 
the xi and X2 fields will be made up of a mixture of many excited states, they 
will have dominant overlap with different states [22,24]. By using the correlation- 
matrix techniques described in Ref. [34] (see also Appendix) , two separate mass 
states are extracted from the xi and X2 interpolating fields. The results from 
the correlation matrix analysis are shown by the filled symbols in Fig. 4, and 
are compared to the standard "naive" fits performed directly on the diagonal 
correlation functions, XiXi a-^id X2X27 indicated by the open symbols. 

The results indicate that indeed the N*{^ ) largely corresponding to the 
X2 field (labeled 'W|") lies above the iV*(i ) which can also be isolated via 
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Euclidean time evolution with the xi Add ("iVj*") alone. The masses of the cor- 
responding positive parity states, associated with the xi and X2 fields (labeled 
'Wi" and "N2" , respectively) are shown for comparison. For reference, we also 
list the experimentally measured values of the low-lying i states. It is interest- 
ing to note that the mass splitting between the positive parity Ni and negative 
parity 7V;f 2 states (roughly 400-500 MeV) is similar to that between the A^*2 
and the positive parity A'^2 state, reminiscent of a constituent quark-harmonic 
oscillator picture. 

Turning to the strange sector, in Fig. 5 we show the masses of the positive 
and negative parity S baryons calculated from the FLIC action [34] compared 
with the physical masses of the known positive and negative parity states. The 
pattern of mass splittings is similar to that found in Fig. 4 for the nucleon. 
Namely, the ^ state associated with the xi field appears consistent with the 
empirical 17(1193) ground state, while the i state associated with the X2 field 

lies significantly above the observed first (Roper-like) ^ excitation, 17* (1660). 
There is also evidence for a mass splitting between the two negative parity states, 
similar to that in the nonstrange sector. 

The spectrum of the strangeness -2 positive and negative parity S hyperons 
is displayed in Fig. 6. Once again, the pattern of calculated masses repeats that 
found for the E and N masses in Figs. 4 and 5, and for the respective coupling 
coefficients [34] . The empirical masses of the physical S* baryons are denoted 
by asterisks. However, for all but the ground state S'(1318), the J^ values are 
not known. 

Finally, we consider the A hyperons. In Figs. 7 and 8 we compare results 
obtained from the A^ and yl"^ interpolating fields, respectively, using the two dif- 
ferent techniques for extracting masses. A direct comparison between the positive 
and negative parity masses for the A^ (open symbols) and yl"^ (filled symbols) 
states extracted from the correlation matrix analysis, is shown in Fig. 9. A sim- 
ilar pattern of mass splittings to that for the N* spectrum of Fig. 8 is observed. 
In particular, the negative parity A^ state (diamonds) lies ^ 400 MeV above the 
positive parity Ai ground state (circles), for both the A^ and A"^ fields. There 
is also clear evidence of a mass splitting between the A^ (diamonds) and A2 
(squares) . 

Using the naive fitting scheme (open symbols in Figs. 7 and 8) misses the 
mass splitting between A^ and yl2 for the "common" interpolating field. Only 
after performing the correlation matrix analysis is it possible to resolve two 
separate mass states, as seen by the filled symbols in Fig. 8. As for the other 
baryons, there is little evidence that the yl2 (triangles) has any significant overlap 
with the first positive parity excited state, yl*(1600) (cf. the Roper resonance, 
A^*(1440), in Fig. 4). 

While it seems plausible that nonanalyticities in a chiral extrapolation [11] 
of Ni and N^ results could eventually lead to agreement with experiment, the 
situation for the yl*(1405) is not as compelling. Whereas a 150 MeV pion-induced 
self energy is required for the A^i, N^ and Ai, 400 MeV is required to approach 
the empirical mass of the yl*(1405). This may not be surprising for the octet 
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fields, as the /1*(1405), being an SU(3) flavour singlet, may not couple strongly to 
an SU(3) octet interpolating fleld. Indeed, there is some evidence of this in Fig. 9. 
This large discrepancy of 400 MeV suggests that relevant physics giving rise to a 
light yl*(1405) may be absent from simulations in the quenched approximation. 
The behaviour of the A\ 2 states may be modified at small values of the quark 
mass through nonlinear effects associated with Goldstone boson loops including 
the strong coupling of the yl*(1405) to Sii and KN channels. While some of this 
coupling will survive in the quenched approximation, generally the couplings are 
modified and suppressed [12,61]. It is also interesting to note that the A\ and 
A2 masses display a similar behaviour to that seen for the H* and S'2 states, 
which are dominated by the heavier strange quark. Alternatively, the study of 
more exotic interpolating fields may indicate the the yl*(1405) does not couple 
strongly to xi or X2- Investigations at lighter quark masses involving quenched 
chiral perturbation theory will assist in resolving these issues. 
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6.4 Spin-| Baryons 



The interpolating field defined in Eq. (40) has overlap with spin-i and spin-| 
states of both parities. After performing appropriate spin projections on the cor- 
relation functions, the masses of the N^ and N-^ states are extracted and 

displayed in Fig. 10 as a function of m^. Earlier results for the N^ states using 
the standard spin-i interpolating field [33,34] from Eq. (16) are also shown with 
open symbols in Fig. 10 for reference. It is encouraging to note the agreement be- 
tween the spin-projected ^ states obtained from the spin-| interpolating field 
in Eq. (40) and the earlier i results from the same gauge field configurations. 
We also observe that the 7V| state has approximately the same mass as the 
spin-projected N^ state which is consistent with the experimentally observed 
masses. The results for the N^ state in Fig. 10 indicate a clear mass split- 
ting between the N^ and N^ states obtained from the spin-| interpolating 
field, with a mass difference around 300 MeV. This is slightly larger than the 
experimentally observed mass difference of 200 MeV. 
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Turning now to the isospin-| sector, results for the Z\7 



and Z\|" 

,3 + 



masses 



are shown in Fig. 11 as a function of m^. The trend of the Z\| data points 
with decreasing m± is clearly towards the Z\(1232), although some nonlinearity 
with rn^ is expected near the chiral limit [11,12]. The mass of the Z\| lies some 
500 MeV above that of its parity partner, although with somewhat larger errors. 
After performing a spin projection to extract the A^ states a discernible, 
but noisy, signal is detected. This indicates that the interpolating field in Eq. (43) 
has only a small overlap with spin-i states. However, with 400 configurations we 
are able to extract a mass for the spin-2 states at early times, shown in Fig. 11. 
Here we see the larger error bars associated with the A^ states. The lowest 
excitation of the ground state, namely the A^ , has a mass ^ 350-400 MeV 



above the A 



3 + 



with the Z\| 



possibly appearing heavier. The A^ state is 



found to lie ~ 100-200 MeV above these, although the signal becomes weak at 
smaller quark masses. This level ordering is consistent with that observed in the 
empirical mass spectrum, which is also shown in the figure. 

The N^ and A^ states will decay to Ntt in S'-wave even in the quenched 
approximation [62]. For all quark masses considered here, with the possible ex- 
ception of the lightest quark, this decay channel is closed for the nucleon. While 
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there may be some spectral strength in the decay mode, we are unable to sepa- 
rate it from the resonant spectral strength. 

The A'' I and A^ states will decay to Nn in P-wave, while A^| and Z\| 
states will decay to Nn in D-wave. Since the decay products of each of these 
states must then have equal and opposite momentum and energy given by 

E^ = m2 + (— 

these states are stable in the present calculations. 



7 Conclusions 

The increasing effort given to the study of the excited baryon and meson spec- 
trum by the lattice community reflects the appreciation that the determination 
of the spectrum provides vital clues to the dynamics of QCD, and the mecha- 
nisms of confinement. The N* spectrum in particular has several features, such 
as the anomalously light iV*(1440) Roper resonance and the J^ = ^ yl(1405), 
that defy a straightforward interpretation within the quark model, and whose 
understanding can help to resolve the competing pictures of hadron structure. 
The impetus to study the N* resonances has also strengthened following the 
observation of the S" = -|-1, 0^ pentaquark state, whose properties remain essen- 
tially unknown, but whose various interpretations offer vastly different pictures 
of baryon structure. 

In these lectures, we have described the techniques required to compute the 
A'^* spectrum in lattice QCD, and given an overview of the current status of 
lattice calculations. The lightest states of both parities for spin 1/2 and spin 
3/2 have been successfully resolved in the quenched approximation to QCD in 
both the nucleon (isospin-i) and A (isospin-|) sectors, and in general the level 
ordering, albeit at relatively large pseudoscalar masses, 111^^ > 500 MeV, is in 
accord with that observed experimentally. At these large pseudoscalar masses 
the spectrum largely follows quark-model expectations. 

It is important to appreciate that most current spectroscopy calculations 
have employed essentially "S-wave" quark propagators. The measurement of a 
wider basis of interpolating operators will be an important element of future 
studies, and the technology to construct such operators with lattice symmetry 
properties has now been developed. An important by-product of such studies 
will be insight into the quark and gluon structure of such hadrons. 

The realisation that hadronic physics at physical values of the pion mass is 
very different from that at 111^^ > 300 MeV, and the consequent need to correctly 
account for the chiral properties of the theory, have been some of the most 
important developments in lattice QCD of recent years. FLIC fermions and the 
development of fermions having an exact analogue of chiral symmetry provide 
the means to attain such pion masses. There are suggestions that such radically 
different behaviour at light quark masses has been seen in the excited nucleon 
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sector. With the development of exact dynamical fat-link fermion algorithms 
[63-65] , FLIC fermions provide tremendous promise for accessing the light-quark 
mass regime of full QCD. 

The continuation to physical values of the light quark masses poses extra 
challenges to the calculation of the resonance spectrum. The excited states are 
no longer stable under the strong interaction at sufficiently light quark masses. 
Even in the quenched approximation, this instability is manifest through non- 
unitary behaviour in the correlators and through additional non-analytic terms 
in the chiral expansion of nucleon masses [13]. The means to study the full spec- 
trum, including the scattering lengths of multi-particle states, is in principle 
known, relying on examining the finite-volume shift in the two-particle spec- 
trum. However, the method is computationally very demanding, requiring the 
measurement of many operators. 

With the realisation of large-scale computing facilities expected over the next 
several years, we can expect many of these endevours to come to fruition. An 
exciting era for baryon spectroscopy on the lattice lays ahead. 
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8 Appendix - Correlation Matrix Analysis 

In this section we outline the correlation matrix formalism for calculations of 
masses, coupling strengths and optimal interpolating fields. After demonstrating 
that the correlation functions are real, we proceed to show how a matrix of such 
correlation functions may be used to isolate states corresponding to different 
masses, and also to give information about the coupling of the operators to each 
of these states (see also Ref. [34] ) . 

8.1 The U + U* method 

A lattice QCD correlation function for the operator xiXj^ where Xi is the i-th 
interpolating field for a particular baryon (e.g. Xj ^^ Section 5.2), can be written 
as 

Q,, = {Q\T{x^Xj)\n) (64) 
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where spinor indices and spatial coordinates are suppressed for ease of notation. 
The fermion and gauge actions can be separated such that S'[C/, f/j, ^] = 6*0 [C/] + 
ipM[U]'il'. Integration over the Grassmann variables ip and ^ then gives 



/X>t/e-^G[^l dct(M[U])H^j[U] 
JVUe-Solu]dct{M[U]) 






where the term Hij stands for the sum of all full contractions of XiXj- The pure 
gauge action So and the fermion matrix M satisfy 

Sg[U] = Sg[U*] , (66) 

and 

CM[U*]C-^ ^M*[U] , (67) 

respectively, where C = C75 . 

Using the result of Eq. (67), one has 

det {M[U*]) = det {M*[U]) , (68) 

and since det(A'/[f7]) is real, this leads to 

det (M[[/*]) = det (M[[/]) . (69) 

Thus, U and U* are configurations of equal weight in the measure 
J 'DUdet{M[U])exp {— Sg[U]) , in which case Gij can be written as 



1 //X>L/e-SG[^ldet(M[t/]){i7y[C/]+i?y[C/*]}' 

2 I JVUe-Solu]iet{M[U]) 



g. = ^ r r-.:::z:z:'r;:^Jr: ''^ '' ] ■ m 



Gf^ ^ tr[r±g,,] , (71) 



Let us define 

where "tr" denotes the spinor trace and I± is the parity-projection operator 
defined in Eq. (15). If tr [rHij[U*]] = tr [r7J*-[t7]], then Cf^ is real. This can 
be shown by first noting that Hij will be products of Dirac 7-matrices, fermion 
propagators, and link-field operators. In a 7-matrix representation which is Her- 
mitian, such as the Sakurai representation, C^^C^^ ~ 7* Fermion propaga- 
tors have the form A'/^\ and recalling that since C M[U*]C^^—M*[U], then 
we have CM-^[U*]C-^ = {M-^[U])* . For products of link-field operators 0[U\ 
contained in _ffjj, the condition 0[f7*] = 0*[C/] is equivalent to the requirement 
that the coefficients of all link-products are real. As long as this requirement 
is enforced, we can then simply proceed by inserting CC~^ inside the trace to 
show that the (spinor-traced) correlation functions Gf, are real. If one chooses 

the Dirac representation, then C^kC~^ — — 7j! and C7oC^^ — 7q. Therefore, 
in the Dirac representation of the 7-matrices, if Hij contains an even number 
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of spatial 7-niatrices with real coefficients, G^ is purely real; otherwise G^ is 
purely imaginary. 

In summary, the interpolating fields considered here are constructed using 
only real coefficients and have no spatial 7-matrices. Therefore, the correlation 
functions Gj are real. This symmetry is explicitly implemented by including 
both U and U* in the ensemble averaging used to construct the lattice correlation 
functions, providing an improved unbiased estimator which is strictly real. This 
is easily implemented at the correlation function level by observing 

M-\{U;}) = [G75M-i({C/^})(G75)-^]* 
for quark propagators. 

8.2 Recovering Masses, Couplings and Optimal Interpolators 

Let us again consider the momentum-space two-point function for t > 0, 

g,j{t,p) = ^e-^P-(/7|x.(t,a^)Xj(0,O)|l2) . (72) 

At the hadronic level, 

X p' ,s B 

X {B,p\s\xj{0,O)\n) , 
where the \B,p' , s) are a complete set of states with momentum p' and spin s 

EEEi^'p''^)(^'p''^i==^- (73) 

p' B s 

We can make use of translational invariance to write 

t;,,(t,p)=^e-P-5]5]^(l2|e^*e-^-X.(0)e^^-^e^^*|s,p',5 

X p' s B 

x{B,p\s\x,{Q)\n) 
= EE^"''"'(^I^«(0)I^'P'*)(^'P'*IX,(0)I^) • (74) 

s B 

It is convenient in the following discussion to label the states which have the 
X interpolating field quantum numbers and which survive the parity projection 
as \Ba) for a = 1, 2, • • • , A^. In general the number of states, iV, in this tower 
of excited states may be infinite, but we will only ever need to consider a finite 
set of the lowest such states here. After selecting zero momentum, p = 0, the 
parity-projected trace of this object is then 

N 

G± (i) = tr[r±g,, (i, 0)] = Y. e"""* A? a; , (75) 

a = l 
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where \f and A.- are coefBcients denoting the couplings of the interpolating fields 
Xi and Xj, respectively, to the state \Ba). If we use identical source and sink 
interpolating fields then it follows from the definition of the coupling strength 
that aJ = {X^y and from Eq. (75) we see that Gfj{t) = [G%{t)]\ i.e., G± 
is a Hcrmitian matrix. If, in addition, we use only real coefficients in the link 
products, then G is a real symmetric matrix. For the correlation matrices 
that we construct, we have real link coefficients but smeared sources and point 
sinks. Consequently, G is a real but non-symmetric matrix. Since G is a real 
matrix for the infinite number of possible choices of interpolating fields with real 
coefficients, then we can take Af and A, to be real coefficients without loss of 
generality. 

Suppose now that we have M creation and annihilation operators, where 
M < N . We can then form an M x M approximation to the full N x N matrix 
G. At this point there are two options for extracting masses. The first is the 
standard method for calculation of effective masses at large t described in Sec- 
tion 5.1. The second option is to extract the masses through a correlation-matrix 
procedure [44]. 

Let us begin by considering the ideal case where we have N interpolating 
fields with the same quantum numbers, but which give rise to N linearly in- 
dependent states when acting on the vacuum. In this case we can construct N 
ideal interpolating source and sink fields which perfectly isolate the N individual 
baryon states \Ba), 

N 



5"-E<^. , (76) 

i=l 

N 

*"=E<"^- (77) 



i=l 



such that 

(B^in^)-'5a;3z"u(a,P,s) , (78) 

(/?|(/)"|B^)=(5„^z"u(a,p,s) , (79) 

where 2" and J" are the coupling strengths of (/)" and (j) to the state \Ba). The 
coefficients uf and ii*" in Eqs. (77) may differ when the source and sink have 
different smearing prescriptions, again indicated by the differentiation between 
z" and 2". For notational convenience for the remainder of this discussion re- 
peated indices i,j,k are to be understood as being summed over. At p = 0, it 
follows that, 

= Afz^e^""*. (80) 
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The only i-dependence in this expression comes from the exponential term, which 
leads to the recurrence relationship 

G%{t)u^^e"^-Gi{t + l)ul, (81) 

which can be rewritten as 

[G^{t + lKlG%{t)u'^^e"^-ul. (82) 

This is recognised as an eigenvalue equation for the matrix [G^{t + l)]^^G*(t) 
with eigenvalues e"'° and eigenvectors u". Hence the natural logarithms of the 
eigenvalues of [G^{t + l)]^^G^(i) are the masses of the N baryons in the tower 
of excited states corresponding to the selected parity and the quantum numbers 
of the X fields. The eigenvectors are the coefficients of the x fields providing the 
ideal linear combination for that state. Note that since here we use only real 
coefficients in our link products, then [G^{t + l)\~^G^{t) is a real matrix and 
so ii" and w" will be real eigenvectors. It also then follows that z" and z" will 
be real. These coefficients are examined in detail in the following section. 

One can also construct the equivalent left-eigenvalue equation to recover the 
V vectors, providing the optimal linear combination of annihilation interpolators, 

«rG^,(i)-e™°<"G±(i + l). (83) 

Recalling Eq. (80), one finds: 

G%{t)u'^ =z"A?e-™°*, (84) 

Gfm =z"re-™°*, (85) 



vr 



i± u\n±i 



«r G%{t)Gt{t) uf = z"z"Af A, e-^"°* . (86) 

The definitions of Eqs. (79) imply 

<" G±(i) Mj" = z^z^e-""*, (87) 

indicating the eigenvectors may be used to construct a correlation function in 
which a single state mass m^is isolated and which can be analysed using the 
methods of Section II. We refer to this as the projected correlation function in 
the following. Combining Eqs. (86) and (87) leads us to the result 

v;- Gk,{t)Ga{t) u'^ a _„^, 

By extracting all N"^ such ratios, we can exactly recover all of the real couplings 
Xf and A, of Xi a-nd %, respectively to the state |-Bq). Note that throughout 
this section no assumptions have been made about the symmetry properties of 
Gj . This is essential due to our use of smeared sources and point sinks. 

In practice we will only have a relatively small number, Af < iV, of interpo- 
lating fields in any given analysis. These M interpolators should be chosen to 
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have good overlap with the lowest M excited states in the tower and we should 
attempt to study the ratios in Eq. (88) at early to intermediate Euclidean times, 
where the contribution of the {N — M) higher mass states will be suppressed but 
where there is still sufficient signal to allow the lowest M states to be seen. This 
procedure will lead to an estimate for the masses of each of the lowest M states 
in the tower of excited states. Of these M predicted masses, the highest will in 
general have the largest systematic error while the lower masses will be most 
reliably determined. Repeating the analysis with varying M and different com- 
binations of interpolating fields will give an objective measure of the reliability 
of the extraction of these masses. 

In our case of a modest 2x2 correlation matrix (M = 2) we take a cautious 
approach to the selection of the eigenvalue analysis time. As already explained, 
we perform the eigenvalue analysis at an early to moderate Euclidean time where 
statistical noise is suppressed and yet contributions from at least the lowest two 
mass states is still present. One must exercise caution in performing the analysis 
at too early a time, as more than the desired M — 2 states may be contributing 
to the 2x2 matrix of correlation functions. 

We begin by projecting a particular parity, and then investigate the effective 
mass plots of the elements of the correlation matrix. Using the covariance-matrix 
based x^/^'df, we identify the time slice at which all correlation functions of the 
correlation matrix are dominated by a single state. In practice, this time slice 
is determined by the correlator providing the lowest-lying effective mass plot. 
The eigenvalue analysis is performed at one time slice earlier, thus ensuring the 
presence of multiple states in the elements of the correlation function matrix, 
minimising statistical uncertainties, and hopefully providing a clear signal for the 
analysis. In this approach minimal new information has been added, providing 
the best opportunity that the 2x2 correlation matrix is indeed dominated by 
2 states. The left and right eigenvectors are determined and used to project 
correlation functions containing a single state from the correlation matrix as 
indicated in Eq. (87). These correlation functions are then subjected to the 
same covariance-matrix based x^/A^df analysis to identify new acceptable fit 
windows for determining the masses of the resonances. 
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